Uniform approximation for the overlap caustic of a quantum state with its translations 
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The semiclassical Wigner function for a Bohr-quantized energy eigenstate is known to have a 
caustic along the corresponding classical closed phase space curve in the case of a single degree 
of freedom. Its Fourier transform, the semiclassical chord function, also has a caustic along the 
conjugate curve defined as the locus of diameters, i.e. the maximal chords of the original curve. If 
the latter is convex, so is its conjugate, resulting in a simple fold caustic. The uniform approximation 
through this caustic, that is here derived, describes the transition undergone by the overlap of the 
state with its translation, from an oscillatory regime for small chords, to evanescent overlaps, rising 
, to a maximum near the caustic. The diameter-caustic for the Wigner function is also treated. 
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CN , I. INTRODUCTION 

>\ _ 

O ■ It is often assumed that the Wigner function VF(x) in the phase plane x ~ ip,q), for a semiclassical WKB 
] state can be approximated by a Dirac (5-function on the corresponding classical closed curve in phase space. This 
^ , simplest approximation does, indeed, lead to reliable expectation values for smooth classical observables. Nonetheless, 
' it is hopelessly inadequate for the description of delicate interference effects that are becoming ever more accessible 
to experiments related to quantum information, either in quantum optics, atom traps, or other quickly developing 
r~~i • technologies. It is then necessary to resort to more refined semiclassical descriptions for the phase space representations 
'"^ of quantum states, such as Berry's uniform approximation for the Wigner function 01 ■ 

A typical interference experiment superposes two modified copies of the same initial state (see e.g. @). For 
instance, in quantum optics, it is easy to achieve the unitary transformation that corresponds to a uniform phase 
. space translation (or displacement). This translated state can then interfere with the original state. In general, the 
^ ' unitary translation operator 



Oh: 



~ exp 



^(lAx) 



exp 



(1) 



acts on the state to produce the new state IV'^) = in strict correspondence to the classical translation: 

■ x 1-^ x-l-^. [l^ Thus, given an arbitrary superposition of a state and its translation, a\tp) + b\'tp^) , with |ap -I- |5p = 1, 
. the probability that this is measured to be in the untranslated state is \a + 6(V'|V'^)P- 

Evidently, measurements of such probabilities (through repeated preparation) supply detailed quantum information 
concerning these initial states. It so happens that the full set of possible overlaps defines the complete phase space 
[ — ' representation, 

o ■ 



X(0=(^(^|T_^I^). (2) 
This is known as the chord function the quantum characteristic function (or the Weyl function as in Q), which is 



X 

' the Fourier transform of the Wigner function: 



The latter can be redefined, following Royer [1,0, as 



(3) 



W{^) ^ (4) 

[Tin) 

where i?x, the Fourier transform of the translation operators, corresponds classically to the phase space reflection 
through the point x, i. e. Xq i— > 2x — Xq. An important consequence is that the phase space correlations 0], for 
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translational interference, 

CiO = Mi^^)\' ^ {27Th)' I d,7c'^^^/S|x(»7)l' = (271^)^1 dxW(x)W(x-0, (5) 

coincide with the autocorrelation of the Wigner function itself. A study of interference phenomena using the quantum 
phase space formalism, besides some properties of phase space distributions can be found in Q . 

In the limit of small displacements, ^ — > 0, the correlations attain their maximal value, C(0) = 1. Increased 
translations reduce them to an oscillatory regime. We shall have be concerned with the correlation of states, whose 
chord function can be semiclassically approximated by 0, Q 

3 3 

where the amplitudes and phases are determined by a classical curve, as will be described in the next section. This 
is similar to the simplest semiclassical approximation for the Wigner function [3| 

Wix) = E a,(x) e'^^^(-)/^ = E W,ix) . (7) 

3 3 

Furthermore, the Fourier relation between this pair of representations is reflected in the reciprocal relation, which 
specifics the centres 

x.(0=J^, (8) 
for the realizations of the vector ^ as a chord of the classical curve, whereas 

determines the chords that have a given centre. Here, 

is the standard symplectic matrix. Typically for an energy eigenstate. the classical curve is a level curve for the 
corresponding classical Hamiltonian, related by Bohr-Sommerfeld [§, [13, El, quantization. In this way ([6]) and ([7]) 
are alternative phase space representation of WKB wave functions. 

For large enough displacements, such that the classical translation of the curve docs not intersect the original 
curve, the phase space correlations are negligible. The transition between this and the previous oscillatory regime 
takes place along a caustic where the chord function attains locally maximal amplitudes and the simple semiclassical 
approximation ([6]) breaks down. The main purpose of this paper is to establish the correct description of this transition 
from the oscillatory regime to the region of negligible overlap through a uniform approximation. It is interesting that 
the caustic region for increased quantum correlations is entirely determined by the geometry of the classical curve 
supporting the quantum state. It will here be assumed that this curve is convex, so that the locus of its diameters, 
i.e. its maximal chords, also defines a closed convex curve. 

The present approximation has some resemblance to the uniform approximation obtained by Berry 0] for the 
Wigner function close to the classical curve. However, the latter is simplified by symmetry constraints that do not 
hold here. Indeed, the present treatment is even closer to the uniform approximation along the caustic of the Wigner 
function far from the curve, which will also be treated. For a start, section 2 reviews the geometrical construction 
of the Wigner function and the chord function for a state that corresponds to a closed quantized curve. In contrast 
to the treatment in [H, the construction of the present uniform approximations cannot be limited to a single WKB 
branch. 

Having defined geometrically the stationary phases for the Wigner function and the chord function, the method of 
Chester, Friedman and Ursell [l3| then supplies uniform approximations for the chord function in terms of the Airy 
function and its derivative in section 3 and for the Wigner function in section 4. In both cases, the asymptotic form 
of these functions for large argument are then connected to the simpler semiclassical forms © and ((T]) . The analysis 
of these uniform approximations close to the caustic furnishes simpler transitional approximations in section 5. 

None of these expressions extends right down to the limit of small chords, because this is another caustic for both 
the Wigner function and the chord function. However, for the eigenstates of the harmonic oscillator (Fock states), 
the approximations for large chords can be compared to a small chord formula specified by a Bessel functional . This 
leads to a discussion in section 6 of the normalization of all these phase space approximations, deriving from the limit 
C(0) = 1. 
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II. CONSTRUCTION OF THE SEMICLASSICAL WIGNER AND CHORD FUNCTIONS 



Before embarking on uniform approximations for the Wigner and chord functions, it is worthwhile to review the 
derivation of the simpler semiclassical formulae presented in the introduction. Thus we also specify explicitly the am- 
plitudes and phases and their geometric interpretation, which are essential ingredients of the uniform approximations. 

The starting point is the generalized WKB expression for the wave function [l^, [ll|, , 



{q\^l)=NY, 



dqdl 



exp 



n 



(11) 



Here we assume that the classical curve is defined as a level curve of the action variable /(x), such that Sj{q,I) 
corresponds to the j'th branch of the generating function for the canonical transformation {p,q) i— > (1,0), (3j is the 
Maslov correction and N is the overall normalization constant. Choosing the arbitrary initial point, go, the action 
branches are defined as 



Sj (g, /) = /" Pj {Q,I)dQ, ^ = 0j, ^ = Pj , 



dq 



so that the amplitude can be rewritten in terms of 

d^Sj{q,I) _ dp, 



dqdl 



81 



81 

dp 



(12) 



(13) 



It will be assumed that the state corresponds classically to a convex closed curve, so that there will always be a 
single pair of branches for the action function. This will hold irrespective of any linear canonical transformation, 
which it may be convenient to make, given that both the chord function and the Wigner functions are covariant with 
respect to such changes of phase space coordinates. It will also be important to recall that the closed curve must 
satisfy the Bohr-Sommcrfeld quantization condition: 



pdq = 27rfi n 



1 



n e Z. 



Expressing the translation and reflection operators within the position representation (see e.g. 
function ^ becomes 



(14) 

I), the chord 



while the Wigner ^ is given by 



W{yi) 



(27rfi) 



(TTfi) 



dq{q^\i^i) {il^i\q )' 



-iip-q/h 



(15) 



(16) 



where, in both equations q^ ~ q± ^5/2. 

In the semiclassical limit, the WKB expression can be inserted, so that we obtain in each case a sum of integrals that 
are dominated by their points of stationary phase. Irrespective of whether these points are sufficiently isolated so as 
to allow for immediate evaluation by the stationary phase method, we need to understand the geometric construction 
that defines them. In both cases, each stationary point defines values of q^ pairs, which are g-coordinates of a pair 
of points, X*, lying on the classical closed curve. In the case of the chord function, each x~ is the intersection of the 
classical curve with its uniform translation by the vector — ^, whereas x+ = x^ + £. This geometry is exhibited in 
Fig. [1] which shows that each chord has two realizations in a convex closed curve. Thus, this construction on a given 
convex curve always specifies a pair of centres xi(^) and X2(^) for each chord ^ that can be fitted in the curve. 

In the case of the Wigner function, instead of a translation, the classical curve is reflected through the reflection 
centre, x. This results in a pair of intersections, x"*" and x~, such that x is the centre for the pair of chords 
^ = ±(x+ — x^). The stationary points are then the q-coordinates for this pair of chords. This geometry is shown in 
Fig. [21 Therefore, there will be at least one pair of chords ±^(x), for each reflection centre in the curve. 

The stationary phase evaluation of each integral for the Wigner function. 



W, 



d^c 



|(p.(<z+),9+)|(P,(g- 



exp 



^[S.iq^ 



,I)-S,iq-,I)-pi, + (3,-P, 



(17) 
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FIG. 1: Geometrical construction for the stationary points of the chord function. The intersections at xi~ and X2~, of the 
closed curve with its translation by — ^, define the pair of realizations of this chord. The g-coordinate of the centres, xi and X2 
of these realizations defines the stationary points. 




FIG. 2: Geometrical construction for the stationary points of the Wigner function. The intersections of the closed curve with 
its reflection through x defines pair of chords centred on x. The g-coordinate of these chords defines the stationary points. The 
area j4(x) determines the phase of the Wigner function. 



can usually be obtained from a single branch of the action function for the closed curve, that is, i = j. The stationary 
phase is half of the area between the curve and its reflection, or the "chord area" as shown in Fig. [2l except for 
a Maslov correction. The only difference between both phases, corresponding to ±^(x), is the sign, so that the 
semiclassical approximation is real 0. 

In the case of the chord function, each stationary phase is given by the construction in Fig. [31 Unlike the Wigner 
function, this depends explicitly on a change of phase space origin, x x + according to the exact formula [7|, 

= (18) 

but it is also covariant with respect to homogeneous linear canonical transformations. 

The amplitudes in the above semiclassical approximations are best expressed in terms of the canonical action 
variable, /(x), that defines the closed curve and is conjugate to the angle variable 9{x) along the curve. If we now 
define the transported action variable, 

/±=/(x±|/2), (19) 
then, generally, the Poisson bracket for this pair of functions. 



(20) 
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(21) 



and it is found that these amplitudes in ([6]) and ([7]) are 

a(x) = |{/+,/-}|-i=«(0. 

The difference between considering the ampHtude as a function on x or ^ depends on which variable is fixed in (j2ip . 
This only changes the sign of the Poisson bracket. Here, the equality of coefhcients for both representations holds for 
the chord and its centre, between a specific pair of points (x^,x+) on the closed curve. 

A neat interpretation for these amplitudes follows from the identification of the action variable /(x) with a classical 
Hamiltonian. Then the closed curve becomes a closed trajectory, tangent to the phase space velocity vector, x, and 



X A x^ 



(22) 



depend on the degree of transversality of the 
, 01 and so they diverge at caustics, where x+ 



as shown in Fig. [D It follows that the amplitudes aj{^) (or aj(x)) 
intersection between the curve and its translation, or its reflection [! 
and x~ are parallel. 

So far, we have presumed that the contribution of each stationary phase point to the integrals for either the Wigner 
function or the chord function, can be obtained by considering a single branch of the WKB wave function. Thus, 
i = j in (|17p and this single branch can always be accessed through a canonical phase space rotation in the simple 
semiclassical approximations above and even for Berry's uniform theory for the caustic that arises in the limit of 
small chords [2|. However, this is not possible in the present treatment of the caustic at maximal chords, where ((2T|) 
also diverges for both representations, because the tips of the stationary chord become turning points for the WKB 
function in this limit. It is thus necessary to study this limit with the aid of phase space coordinates, such that it is 
the cross term between the pair of different branches, of the action function, which have stationary phases in (|17p . 
For this reason it will be important to take care of the phase relation between the branches across a turning point: 



{q\iJi)=N 



d^S+{q,I) 



dqdl 



exp 



J^S+{q,I) 



in/2 



d^S^iqJ) 



dqdl 



exp 



(23) 



where 7r/2 is the Maslov phase and N is the normalization constant. Inserting the above wave function in (|15p . we 
obtain four integrals to evaluate. An appropiate choice of the coordinate axes orientation, in which the both of chord 
realizations are crossed, reduces the chord function for (|23p to the single integral 



x(0 



2TTh 



dq 



a/ 

dp 



3T 

{p+{q+),q+) — {p_{q-),q~) 



dp 



Xexp(-[S+{q+,I)-S-iq-,I)-ip-q] +^ 



(24) 



Here q^ = g±^^/2. As stated previously, this geometry can be guaranteed by a phase space rotation. The evaluation 
of the chord areas in the Wigner function and the chord function for this geometry is discussed in Appendix A. 



III. UNIFORM APPROXIMATION FOR THE CHORD FUNCTION 



Let us allow the pair of these stationary points of ([M)) to coalesce for the chord that corresponds to a diameter 
of the closed curve, i.e. a maximal chord, at which the semiclassical amplitude (pij) diverges. In the present case of 




FIG. 3: Several geometrical interpretations for the phase for the semiclassical chord function in the simpler approximation for 
a WKB function, considering one branch. The stationary phase itself is determined by the shaded area in the three casesQ- 



6 




FIG. 4: Geometrical interpretation for the semiclassical amplitudes of both the Wigner function and the chord function. 




FIG. 5: The difference in chord areas, A12, coincides with the area between the closed curve and its translation. This area 
shrinks to zero when { becomes a diameter and its conjugate chord = xi — X2 — > 0. The simplest choice of phase space origin 
is at midpoint between xi and X2. 



a convex curve, these diameters are the locus of a fold caustic with no higher singularity. This is simpler than the 
geometry for the corresponding caustic of the Wigner function, studied in the following section. 

We also simplify the calculation by an appropriate choice of origin, in view of the simple translation property (jlSp 
of the chord function. This ideal origin lies midway between the centres for the pair of chord realizations, shown in 
Fig. m The pair of stationary points qi and (72 of the integrand in ([24|l are solutions of the equation: 

p+{q + ^g/2)-p^{q-^g/2) = ^p, (25) 

which are identified as the position coordinates of xi and X2 in Fig. [S] According to Appendix A, the Bohr-Sommcrfcld 
quantization rule p4)) leads to the corresponding phases as 

5+(gi+C,/2)-5-(<Zi-e,/2)-epgi = A(0+xi A^, (26) 
S+{q2 + eg/2) - S-{q2 - (,/2) - ^^92 = ^2(C) + X2 A ^ (27) 



Instead of evaluating (|24|) by stationary phase, this integral is mapped onto the standard form for a fold diffraction 
catastrophe [Tst : 



X{0 = 



2nih 



exp 



where we have defined 



T.A = Ai +A2, 



iTA 
~2h 



dug{u;i) e 



A12 _ ^1 - ^2 + r? A $ 
2ft ~ 2h 



(28) 



(29) 



The action difference, Ai2^ is the main ingredient in the present application of the method of uniform approximation 
(13I . [TBj . Its geometric definition is the area between the closed curve and its translation, as shown in Fig. O At the 
caustic, the chord ^ = is maximal and its conjugate chord Q: = Xi — X2 — > 0. 

The above integral would define an Airy function [l7| , Ai(C), if g{u](^) = 1, but here the mapping between the 
variables q^ u, respectively in (|24|) and (pS)) . leads to 



^i9i^)+U2)^i9iu)-U2) 



dq 
du 



(30) 
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The approximation now consists in replacing this by a Hnear function which coincides with g{u] ^) at the stationary 

xiq = qi and q{—C^) = X2g = 92 • The Jacobian of the mapping q ^ u 

2hC^ 



points, ±^2 of (pS)) . These map onto ^(C^' 
at the stationary points is specified by 



so that 




-e,/2)-|(g2-W2) 



2n0 



Thus, recahing the definition of the transported action (fTO)) for each chord reahzation, we now define 



"-|{/2+,/2~}| " and 



so as to obtain 



12 



C3 



-A/ 



12 



(31) 
(32) 

(33) 

(34) 
(35) 



as an approximation to (|30p . which is hnear with respect to u and has the correct values at the stationary points. 
Thus, the integral for the chord function becomes 



ih V2 
7V2 1 



exp 



exp 



du C*S/ 



12 



-A/ 



12 



C^S/i2Ai(-C)--^A/i2Ai'(-C), 



(36) 



where Ai' is the derivative of the Airy function. Finally, recalling the definition of the intermediate variable, C, in 
terms of areas (j29p . we obtain the full unitary approximation for the chord function: 



X(C) = exp 



[IA2] 



12 



-Ai 



Ah 



All 



AI 



12 



-Ai' 



3^12 



(37) 



The Airy function, Ai(— ^), oscillates with increasing amplitude as its argument increases and then decays expo- 
nentially for positive values of — C. The maximum amplitude, just bellow the origin, indicates the singularity of the 
simple semiclassical amplitude at the caustic. In this region the second term depending on Ai' (— C) can be neglected, 
but it is necessary to obtain the correct limiting behaviour in the oscillatory region where the simple semiclassical 
description is valid. This term is a new feature in comparison with the uniform approximation for the Wigner function 
for small chords 0], where it is absent because of the reflection symmetry. The other novel feature is the oscillatory 
phase proportional to TiA along the caustic, which will be important to separate the contribution of each realization 
at the oscillatory regime. Indeed in the case when the chord areas are significantly greater than Planck's constant, 
the functions in p7p can be replaced by the asymptotic forms for large negative values (l7j , 



Ai(-a;) 



■ cos 



2 3 TT 



and Ai'(— x) 



: sm 



2 3 TT 

t:X^ 



in order to obtain the correct form in the oscillatory regime as in ([6]): 



X(C) 



exp 



exp I 



\{it.ii}\^ 



IU2+,^2"}I^ 



(38) 



(39) 



This result corresponds to the sum of contributions for each chord realization in the simpler stationary phase approx- 
imation of the semiclassical chord function in where it was considered that each chord realization lies on a single 
of the WKB wave function. 
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FIG. 6: Wigner function catastrophes for a closed curve. The curve itself is a fold {short-chord catastrophe). Another singular 
curve lying in the interior is composed of folds and cusps {long-chord catastrophe). The catastrophe points correspond to centres 
of the diameters of the classical curve. 



IV. UNIFORM APPROXIMATION FOR THE WIGNER FUNCTION 



As the Wigner function is evaluated at a point x that lies further and further inside a convex closed classical curve, 
a caustic will be crossed. This is typically a cusped triangle^], in which there are three chords for each centre, as 
shown in Fig. [SI Generic paths to its interior enter the triangle through a fold caustic joining two cusps. Such a 
caustic point {x^} is the centre of a diameter, as show in Fig. [SI as well being the centre of the other chord, which 
was followed in the continuous path through the caustic. An example of the full fringe pattern where the regions 
characterized by one or three chords are clearly discernible can be found in [Sj . 

The uniform approximation to be derived here concerns the pair of chords, and ^2 i^ Fig- that are born from 
the diameter at the caustic. Away from the cusps, the contribution of the separate chord, ^3 in Fig. [3 can still be 
evaluated by stationary phase. This is just a simple scmiclassical contribution, 

W,{k) = ^ |{/(x + $3/2), /(x - ^3/2)}r ^ cos - j) , (40) 

which could be obtained from a single WKB branch, i ~ j in (I17p . This could also be derived from a cross- branch 
i ^ j (by rotating the phase space coordinates) from the same integral that furnishes the joint contribution of the 
pair of chords, and which coalesce at the diameter 

This crossed chord picture is essential for the uniform approximation. As in the theory for the chord function, we 
now map the integral in the region corresponding to the chords and ^2 onto the diffraction integral, as in (|28p . 

iy(x) = 25Re^exp f z— j y dz g{z;^) e'l- '^''\ + W^i^). (41) 

Here the parameters of the transformation are given by 

E^^Ar(x) + A2(x), ^Ci(x) = ^ = ^iW^, (42) 



in which A12 is the symplectic area bounded by the closed curve and its reflection between the ends of and ^2? 



show Fig. [T] The procedure is straightforward as for the chord function, so that recalling that /,*(x) = /(x ± £,j/2), 



we now define 

A/i2^i{/i+,/r}r^-i{/2+,^2"}r^ and E/12^ |{/+,/r}r^ + i{4,^2"}r^- (43) 

Hence, the linear approximation for the amplitude in (|4ip is 

(44) 



5(z,x) = V2fi (^C^E/12 + ^^^12) > 



so that; the Wigner function is given by 



FIG. 7: Stationary chords near to long-chord catastrophe of the Wigner function, displaying the compositions of the area Aj 
for each stationary chord. The phase difference for coalescent and ^2 chords is the symplectic area A12 limited by their tips. 



W{x) = 2V2N^ 



S/i2sin(f| 



fit 



Ai - 



3A 



12 



4fi 



Ai' 



3Ai2 

Ah 



M^3(x). 



(45) 



Again, as for chord function, the behaviour near the caustic is described correctly in terms of the Airy function and 
its derivative. The fact that the Wigner function also has cusps, which are of higher order than the catastrophes of 
the chord function, adds an additional term. Crossing the fold caustic, the coalescent chords disappear and only the 
additional term remains, coinciding with the simpler stationary phase approximation far form the causticQ, although 
the normalization constant must be re-evaluatcd. 

For regions where A12 S> fi, we can replace the Airy function and its derivative by their asymptotic forms for large 
negative values, obtaining 



W(x) 



i{/+/r}l^ 



(46) 



which is a sum of oscillatory terms, each one with a different phase, as in ([7]). This asymptotic form is a superposition 
of the individual stationary phase approximations for each stationary chord. 



V. APPROXIMATIONS FOR THE TRANSITIONAL REGIONS 



The uniform approximation (j37p is not explicitly resolved very close to the caustic, because A12 — s- as the caustic is 
approached, whereas E/^^ 00. The classical curve can be approximated by parabolae in both the neighborhoods of 
the tips of the realization of This equates the amplitude associated for each stationary point, i.e. S/^^ = 2{/+, /"} 
and A/^^ —> 0. Thus, the term of the derivative of the Airy function in ((37|) cancels near to the caustic. 

To obtain an explicit expression for the transitional chord function, we start by recalling that the action variable, 
/(x), can be interpreted as a Hamiltonian, such that the classical curve is a trajectory, i.e. the level curve /(x) = X. 
Considering x as the centre of a chord rj that conects two points of the curve, x^ and x;,, we obtain as a first 
approximation, x^ ~ x^ + rxa, if x is very close to this curve. Then the action can be expanded as (see Appendix B 
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FIG. 8: The area Si for a point very close to the level curve, /(x) = T. The tips of the chord r] evolve by the action of the 
Hamiltonian /(x). As a first approximation, this evolution is linear. 



inRcf. [|) 



(47) 



where 3x is the Hessian matrix of / at the point x. This quadratic Hamiltonian generates the assumed linear motion. 
On the other hand, the area between tj and the curve (Fig. [8|) is given by [1]. 



S'l(x) ~ — T^X^xX. 



(48) 



Thus, noticing that the centres of the reaUzations of r) are = ±^/2, the symplectic area A12 in Fig. [5] may be 
obtained as 



(49) 



where t± is the time of flight between the tips of each realization of 77, under the action of the Hamiltonian /(x). 

Recalling that Xi and X2 are the midpoints of the realizations of ^, the Poisson brackets of the action in the 
amplitudes in ()37|) will be given by the symplectic products. 



x+ A xi = {/+, /i } ~ -{I+, I^] = X2- A x+ 



So that defining the 'accelerations' 0, [l^, xi and X2, as 



d_ 
dx 



such that x^ ~ + t^x-i and Xj^ ~ X2 + t_X2 , we obtain 



{It, in 



The times r_|. and t_ can be found using (|T7)) . so that becomes 



[X+3iX+]i [X-3 iX-]i 
Thus, we have all the necessary ingredients to obtain the transitional form of the chord function. 



X(l) 



TfUi [T+X-3iX+ +T_X+3_|X-]^ 



xAi 



-23 



[X+JiX+]3 [X-3_iX-]^ 



(50) 



(51) 



(52) 



(53) 



(54) 



If the curve has a local symmctrj^ of reflection with respect to the origin, the hessian matrices will be equal, i.e. 
3g = 3_( and the velocity vectors X+ = — X^. Here we recall that the origin depends of the chord ^, since it has 
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been chosen to be the midpoint between the centres of its reaUzations on the closed curve. Thus, the transitional 
chord function reduces in this simple case to 



(27r)i/ii 



[X+J^X^ 



(55) 



For a diameter the argument of the Airy function cancels, because /(±^^/2) = I, and TiA/2 is the chord area of 
Thus, the transitional approximation remains finite, and at the caustic it is close to a local amplitude maximum. 
We can follow a similar procedure for the Wigner function. Defining |, as the average between the stationary chords 
and together with the pair of phase space points, = x± ^/2, we obtain 



4 sin 



[Tly+3y+y+ - T^y 3y-y ]i 



TT^hs [r^y+ly-y - r+y 3y+y^ 



xAi 



-23 



[y^'3y-y-]^ 



(56) 



Unlike the chord function, this transitional form is undetermined at the caustic, in the particular case that the state is 
invariant with respect to a reflection symmetry, because then the entire caustic colapses to a point, where the Wigner 
function has a catastrophe of higher order. 

VI. LONG-CHORD REGIME FOR FOCK STATES 

Now we consider the excited states, |n), of the one dimensional harmonic oscillator whose classical manifold is a 
circumference centred at the origin, i.e. it has reflection symmetry. The quantization condition (|14|) for these circles 
defines 



7r(p^ +g^) = 2TTh ( n 



The exact chord function is given by 



(57) 



(58) 



where i„ is a Laguerre polynomial. For small chords, ^ ^ fi. 



n 



(59) 



gives a good approximation for the chord function [7|, where Jq, is the Bessel function of order zero. Then, the 
asymptotic form of the Bessel function for large values [l3| leads to 



(60) 



We can compare this result with the oscillatory regime ([55)1 . First, due to symmetry, the area of one chord realization 
is complementary to the area of the other and a simple integral gives the semiclassical phase as 



^i(0 = 2^X-|eK/2X- ( i| 



4X arcsin 



iei/2 
^/2f 



(61) 



Moreover, symmetry equates the Poisson bracket for each chord realization, so that terms containing the derivative 
of Airy function will cancel. The Poisson brackets are then evaluated as 



U+,/-} = ICK/2i-(f ) , 



(62) 



12 



so that, the uniform approximation for the chord function becomes 



2X- 1^ 


"■<- 


'3 Ax' 


') 


4 




m 





It follows that the asymptotic behaviour of Airy function, extrapolated for small chords is 



2xiei 



2fi 



(63) 



(64) 



Note that, to lowest order, the argument in the Bessel function ([55)) . v2X|^|, is one half of the complementary area 
to the intersection between the circle and its translation. 

Thus, the asymptotic limit of the chord function for the long-chord caustic reproduces the chord function of small 
chords for Fock states, in an intermediary region. Furthermore, we immediately obtain the normalization constant as 



2n 



(65) 



This is an alternative derivation to Berry's 0- We can replace this value in ([46]) outside the caustic, where the two 
coalescent chords disappear, so as to recover the simpler stationary phase approximation for the Wigner function 



VII. DISCUSSION 



We have shown that the behaviour of both the Wigner function and the chord function near a maximal chord 
singularity can be described by the Airy function and its derivative. Although, the latter becomes negligible for 
points very close to the caustic, it adds an important contribution to the expansion of the semiclassical distributions 
in the oscillatory region, coinciding there with the simpler stationary phase method. The shape of the diameter-caustic 
is different in the phase space of centres, x, where the Wigner function is defined, and in the space of chords, ^. In 
the latter case, the caustic is located on the locus of diameters, ^jj, maximal chords of the original closed complex 
curve. This diameter caustic is symmetrical with respect to the chord origin and it is also convex. If the assumption 
of convexity is relaxed, the symmetry will be preserved, because is also a diameter, but the simple fold caustic 
may then exhibit higher singularities. This is the case for the diameter-caustic viewed in the phase space of centres 
xd = x(^£i). This caustic of the Wigner function has cusps even in the case of a convex quantized curve 01 • 

The pair of caustics that concern us may be viewed as alternative projection singularities of a single (lagrangian) 
surface in a doubled phase space, X = x x+, i.e the product torus for the pair of quantized curves /(x*) = T. 
The alternative description of double phase space in terms of the centres, x = (x+ -|- x~)/2 and ^ = — ^~ leads 
to a description of the double torus that no longer factors Q. The points X^i = {x{^d)t^d) oh the double torus 
project singularly onto both these double phase space coordinate planes. Thus the two-dimensional tangent plane to 
the torus is completely determined by the pair of tangent vectors: (x(^^),0) and (0,^^)). [l^ 

Let us summarize the behaviour of the chord function for any translation: a) a maximum at the origin; b) an 
oscillatory regime, obtained as a superposition of stationary phase terms for each chord realization; c) a region near 
to the maximal chord, i.e. the diameter expressed in terms of the Airy function and its derivative, where the 
amplitude is again maximal and finally d) an evanescent region for chords longer than diameters (also described by 
the Airy functions). 

We have considered only pure states, for which the phase space correlation ([5]) is given by the square modulus 
of the chord function. Thus, using the asymptotic form of the uniform approximation for the chord function at 
the semiclassical regime, we find that, in the oscillatory region far from the caustic, the phase space correlation is 
approximately 



^4 



{/+,/r}"^ + {/+,/2"}'' + 2{/i*-,/r}~M4,^2"r^sin(:^ 



(66) 



i.e. a pair of classical terms associated to each chord realization and a term that represents their interference. 
This formula corrects the semiclassical phase space correlation presented in Q, which also provides a semiclassical 
interpretation for the invariance of the correlation with respect to Fourier transformation. The important point is that 
the maximal reach of the phase space correlations (also the correlations of the Wigner function) in the neighbourhood 
of a diameter of the quantized curve is just = |x(^)P in the simple form given by (|55p . 
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APPENDIX A: CROSSED CHORDS 



Redefining Si = S+ and Sj = S- in (fT7|) , the phase in the integrand for the Wigner function evaluated at the point 
X = {p, q) is 

fiF(Q; x) =S+{q + Q/2, /) -S^iq- Q/2) - pQ. (Al) 

Let go E^nd qi be the turning points on the closed curve, with qo < qi. Choosing the g-axis to pass through qo, 
according to Fig. [51 then at the stationary phase points, Q = ±^g, we have 

S+iq + ^q/2) = / p+(Q)dQ = ai, (A2) 

rqi /•'3-5g/2 

S'_(g-Cg/2) / p+{Q)dQ+ p^{Q)dQ ^ ai + a2 + as + [-as + a4]. (A3) 

Jqo Jqi 

The stationary chord, ^, is here assumed to be crossed, i.e to have both its tips on different action branches. Then, 
the phase becomes 



hF{^q;Ti) = ai - [ai + 02 + a^] - p^g = -[02 + 04] - p£^q ^ a - j) pdq. (A4) 

Defining ^(x) as the area between the stationary chord and the closed curve, leads to 

HFi^g;^) = A(x) ~ j>pdq = A{^) - 21T (^i + ^ . (A5) 

On the other hand, the phase of the integrand in ([M)) that defines the chord function for the same geometry is 

hF^{Q, = S+{Q + ij2) -S-iQ- ij2) ~ ^^Q. (A6) 

The stationary points are position coordinates g, but their geometrical interpretation provides the respective momen- 
tum coordinates, p (as discussed in sec. 2). Again we denote x = {p,q), thus the stationary phase depends of x 
instead of only g. Choosing the g-axis in the same way as before, we obtain 

fiFx(x,0 ^S+{q + Cg/2)-S.{q-^,/2)-p^,+p^,-Cpq, (A7) 
= S+{q + ^J2) -S-{q- ij2) - + X A e (A8) 
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Defining A^iO the symplectic area between the closed curve and the reaHzation of ^, centred on x, as in Fig. [9l 
leads to 

fiFx(x,e)=^x(e)+xAe-27r(^n+i^. (A9) 

In order to implement the uniform approximations in sec. 3 and 4, we can ignore the additional term f pdq in (|A5p 
and (|A9[) . because it is quantized. It should be recalled that the quantization of the closed curve implies that the 
areas, A(x) and A^iO their respective complementary areas, ^'(x) e A'j^{£,), satisfy the rule 



A(x)+A'(x) =^x(0+^x(C) ^ j pdq = 2Tihl^i+^ , nGZ. 



(AlO) 
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